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ABSTRACT 

We consider the dynamical evolution of small-scale substructure in clusters within two extreme alternate 
scenarios for their possible origin: 1) the accretion of groups (or small clusters) on quasi-radial orbits, and 
2) the merger of clusters of similar masses, followed by the decoupling of their dense cores. Using simple 
analytical arguments and checking with numerical computations, we show that objects are destroyed by 
the tidal field of the global cluster potential if their mean density is small compared to the mean cluster 
density within the radius of closest approach of the group or detached core. Accreted groups and small 
clusters are thus tidally disrupted in one cluster crossing. Since the cores of clusters are much denser than 
groups, they are considerably more robust to tides, but the least massive are destroyed or severely stripped 
by tides, while the others are brought to the cluster center by dynamical friction (and subsequently merge) 
in less than one orbit. The longest lived substructures are detached cores, roughly ten times less massive 
than the cluster, starting in near-circular orbits beyond 1 h~ x Mpc from the cluster center. 
Subject headings: galaxies: clustering — galaxies: kinematics and dynamics — cosmology: theory 



1. INTRODUCTION 

i Clusters of galaxies are often complex structures h ar- 
ising two different kinds of probably related substruc- 
^re. First, large-scale substructure, appearing as a bi- 
yioda.l grand design, has been reported in roughly 30% of 
flusters from the study of projected galaxy positions, red- 
snifts, and X-ray emission maps (e.g. Dressier k Sheet- 
man 1988b; Jones k Forman 1992; Escalera et al. 1994). 
^Second, small-scale substructure, in the form of density en- 
hancements around luminous galaxies seems also present 
the central regions of some otherwise apparently re- 
ixed clusters. Such small-scale substructure is appar- 
•it in the Coma cluster (Fitchett k Webster 1987; Mel- 
fer et al. 1988; White, Briel, k Henry 1993), and of- 
«^u|n inferred from gravitational optics analyses of giant 
CTfcs (Miralda-Escude 1993; Kneib et al. 1993). A frac- 
(Thq n of ~ 50% of large-scale monomodal clusters could 

f bor small-scale substructure according to the statisti- 
study by Salvador-Sole, Sanroma, k Gonzalez-Casado 
93, hereafter SSG) of Dressier k Shectman's (1988a) 
iple of clusters, confirmed by the analysis of correlation 
ong galaxy positions within those clusters (Salvador- 
Sole, Gonzalez-Casado, k Solanes 1993, hereafter SGS), 
and from wavelet analysis (Escalera et al. 1994). 

The presence of large-scale bimodal substructure 
strongly suggests that clusters are caught in the process 
of merging. Thus, the amount of clusters harboring large- 
scale substructure is closely related with the rate of for- 
mation of massive cosmic structures. Following this idea, 
three recent studies (Richstone, Loeb, k Turner 1992; 
Lacey k Cole 1993; Kauffmann k White 1993) have used 
the observed fraction of large-scale substructure in clus- 
ters to constrain the value of Slo (basically guessing the 
survival times of substructures), and while the first two 
concluded that Slo > 0.5, Kauffmann k White (1993) re- 
marked that the survival time of substructure is too un- 
certain to conclude on Slo- In this Letter, the dynamical 
survival of observed small-scale substructure in clusters is 



assessed with attention to its possible origin. 

In the framework of hierarchical clustering, small-scale 
substructure could appear in apparently relaxed clusters in 
two ways. First, this substructure could simply be a group 
or small cluster accreting onto the cluster. A more sub- 
tle idea (Gonzalez-Casado, Solanes k Salvador-Sole 1993, 
hereafter GSS) is to have two clusters of similar mass 
merge together, and in that process, their central cores 
detach from their main body. This type of dynamical evo- 
lution is seen in simulations of galaxies merging together: 
their extended dark halos merge before their more bound 
visible components, (e.g. Barnes 1992). Indeed, the vi- 
olent relaxation induced by the merger generates rapid 
variations in the global potential, subjecting different la- 
grangian elements to rapidly varying tides. Although a 
detailed study of tidal processes during merging is well be- 
yond the scope of this Letter, it is safe to assume that, if 
the initial clusters have cuspy cores, these are indeed likely 
to survive the tidal processes during the violent relaxation 
phase occurring during the merger. Both accretion and 
merging scenarios are studied below. 

2. DYNAMICS OF IMMERSED SUBSTRUCTURE 

Groups or small clusters accreting onto larger clusters 
as well as those detached cores that follow elongated or- 
bits within a cluster (both subsystems are hereafter called 
clumps) are subject to a tidal shock from the global po- 
tential of the cluster as they pass by its center. 

The tidal acceleration of a galaxy in the clump can be 
written as 
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where R and S are the positions of the clump center and 
one galaxy within the clump in the frame of the cluster 
center, respectively, R p is the orbit pericenter, and where 
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with f(R) = M(R)/M(R p ), while S and R are unit vec- 
tors along S and R, respectively. Integrating over time, 
assuming that clumps remain stationary, and working in 
cartesian coordinates, one can obtain the velocity incre- 
ment of galaxies in a clump. Then, averaging these in- 
crements over the galaxy distribution of the clump and 
assuming an impulsive approximation (Spitzer 1958; see 
also Mamon 1987) leads to a gain in internal energy 
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where m s and r s are the mass and radius of the clump, 
respectively, V p is the clump's velocity at pericenter, and 
Aw is a dimensionless velocity increment: 
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with T the time spent by the clump between two consec- 
utive apocenters. Writing the clump's internal energy as 
U = —jGm 2 /r s , one finally obtains 
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where p c \ and p s are the respective mean densities of the 
cluster and substructure, while V c i TC is the circular veloc- 
ity of the cluster. Hence, the relative internal energy gam 
scales as the ratio of the mean densities of cluster within 
pericenter to substructure (see also Mamon 1992). Equa- 
tion (3) implies that groups or small clusters that penetrate 
deeply into the cluster will suffer tides more severely than 
detached cuspy cores, since the mean density of a cluster 
core is much higher than the mean density of a whole group 
or cluster. Only groups as compact as Hickson's (1982) 
compact groups may be robust to tidal destruction, since 
they appear as dense as the cores of rich clusters. 

We integrate numerically the orbit of the clump in the 
cluster potential, including the loss of orbital energy of 
the clump by dynamical friction against the dark matter 
particles of the cluster. The force of dynamical friction is 
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where xo = 2~ 1 l 2 v / a c \, p c \(R) and <r c i(i?) are the mass 
density and ID velocity dispersion of the cluster at a dis- 
tance R from its center, and A = Pmax/Pmin, is the ratio 
between the maximum and minimum impact parameters 
for which collisions between the subsystem and the light 
background particles are efficient. Here, we assume either 
1) Pmax = R and p mm = Gm s /v 2 for clusters and clumps 
with a singular isothermal density profile (Tremaine 1976; 
Lacey & Cole 1993); or 2) p max = max(i? c , R) and p mm = 
r c for clusters and clumps with a flat core of radius equal 
to R c and r c , respectively (Tremaine, Ostriker & Spitzer 
1975). 

The tidal acceleration (eq. [1]) is evaluated along the or- 
bit, from first to second apocenter. Velocity increments are 
calculated by numerical integration of equation (2) and the 



relative internal energy gain is computed through equation 
(3). A subsystem with |A?7/i7| > 1 after one orbit is con- 
sidered to be disrupted within a semi-orbital time, T/2. 
Indeed, clumps destroyed by tidal forces will dissolve in 
a timescale shorter than their internal semi-orbital time, 
which itself is shorter than T/2. Clumps surviving cluster 
tides and very massive clumps will have their orbits decay 
by dynamical friction towards the cluster center. In these 
cases, the survival time is set equal to the orbital decay 
time, defined as the time for the clump to definitely settle 
within the cluster core radius. 

3. ACCRETED GROUPS 

From spherical cosmological infall, a group accreted on 
a radial orbit by a cluster and reaching its pericenter at 
present time, to, had an apocenter that was the clus- 
ter turnaround radius at time to/2 (Gunn & Gott 1972): 

1/3 

R a = (2GM c \t 2 /% 2 ) , where M c \ is the present clus- 
ter mass within the radius, R\ = R a /2, of a circular or- 
bit with period equal to to- This roughly approximates 
the size of the region which may be relaxed in the clus- 
ter (White 1992). The cluster crossing time is defined as 
t CT = R } /V cm (Ri) =V(2tt). 

In this accretion scenario, the cluster is assumed to have 
a cuspy density profile, p c \ ~ R~ 2 , consistent with grav- 
itational optics analyses (e.g. Hammer 1991). We adopt 
the same cuspy density profile for infalling groups or small 
clusters, and their radii are inferred in the same way as the 
cluster radius. This contradicts the reported absence of 
central concentration in the density profile of small groups 
(Walke & Mamon 1989), therefore, the tidal disruption 
times derived for groups will be upper limits, since we are 
overestimating their central concentration and, therefore, 
their ability to survive tides. 
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FIG. 1. Relative internal energy increment after a pericen- 
tric passage of accreted groups of internal velocity dispersion <r s , 
with apocenter R a = 2.2 h — 1 Mpc. The initial apocentric velocities 
are 250kms- 1 (s quarts) and 50kms 1 [triangles), corresponding 
to specific angular momenta S = (J / J c irc[E]) 2 — 0-2 and 0.01, re- 
spectively. The mass fraction, /i = m s /M c \, is given in the top axis. 
The horizontal dashed- dotted line is for critical disruption. 
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Figure 1 shows the relative internal energy gains of ac- 
creted groups and small clusters of different velocity dis- 
persions (or mass m s = fiM c \) and initial circularity pa- 
rameters (defined as the squared ratio of the initial an- 
gular momentum to the angular momentum of the cir- 
cular orbit of same energy, Merritt 1985). We adopt a 
ID cluster velocity dispersion of u c \ = 760 km s -1 , yield- 
ing Ri = 2- 1 / 2 7r- 1 (7 c ito = 1.1 h' 1 Mpc and M cl = 
2 1 / 2 (7rG)- 1 (T^ = 3.1 x I0 14 h- 1 M Q , for fi = 1. Fig- 
ure 1 clearly shows that infalling groups and small clus- 
ters are tidally disrupted at first passage through the core 
of a rich cluster, regardless of their initial orbital parame- 
ters, since they gain considerable more energy than their 
binding energy. This result is independent of the values 
assumed for Cto and u c \ > 400 km s -1 . A similar result has 
recently been reported (Capelato & Mazure 1995) through 
use of TV-body simulations (with a non-cuspy cluster po- 
tential). Groups accreted on nearly radial orbits are more 
easily disrupted because their pericenters are closer to the 
cluster center. 

4. DETACHED CORES OF MERGED CLUSTERS 

We now evaluate numerically the survival time of a 
clump, initially immersed in a cluster, similarly as for ac- 
creted groups, except that: 1) the cluster density profile 
is flat at the center, since cores have detached, and 2) the 
initial distance to the cluster center and the initial specific 
orbital angular momentum are free parameters, since, af- 
ter violent relaxation, the clumps are left orbiting with 
random position and velocity direction. Since violent re- 
laxation leaves a phase space distribution independent of 
mass (Lynden-Bell 1967), the typical initial velocity V% of 
clumps will be the 3D velocity dispersion of galaxies in the 
cluster. 

We adopt, for the cluster, a Modified Hubble number 
density profile n(R) = «o/[l + (R/Rc) 2 ] 3 / 2 and constant 
galaxy velocity dispersion a g , both consistent with obser- 
vations. Assuming isotropic hydrostatic equilibrium, the 
gravitational potential of the cluster mass distribution is 
then $(i?) = (3/2)cr 2 ln[l + (R/R c ) 2 ]. The resulting (es- 
sentially dark matter) background mass density and ID 
velocity dispersion are 
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where X = R/R c , R c = 0.25 Mpc is the cluster core 
radius, and where <r c \ is obtained by solving the isotropic 
hydrostatic equation for the background mass. The same 
model has been assumed for the detached core, but with a 
smaller core radius, r c = 0.06 Mpc, following the op- 
tical properties of clumps obtained by GSS from a simple 
clustering model that reproduces the statistical properties 
of the observed small-scale substructure in clusters. Ac- 
cording to SGS and GSS, the small-scale substructure in 
clusters occurs in clumps smaller than 0.3/i _1 Mpc, and 
seems to harbor typically about 10% of the galaxies in the 
central l/i -1 Mpc region of clusters. 
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FIG. 2. Relative internal energy increment of detached cores 

of size r s = 0.2 h~ 1 Mpc vs. initial distance to the cluster center, 
R l: for different mass fractions /i = m s /M c \. Triangles represent 
radial orbits, and squares are for orbits with specific angular mo- 
mentum S = 0.25. The horizontal dashed- dotted line is for critical 
tidal disruption. Only cases with at least one secondary apocenter 
are shown. 

Figure 2 shows the relative gain in internal energy af- 
ter a pericenter passage for detached cores on elongated 
orbits, with mass m s = fiM c \, for r s = 0.2 /i -1 Mpc (we 
have tried different values of r s ). Massive clumps with 
fi > 0.05 and r s < 0.25 Mpc survive the cluster tidal 
field, whereas low-mass clumps do not. For nearly circu- 
lar orbits, the tide is not impulsive but roughly station- 
ary, and the clump may be tidally limited to some radius. 
Severely limited clumps cannot be observed as small-scale 
substructure. The lower limit of the tidal radius is ob- 
tained when the clump's orbit is circular, and its internal 
rotation is phase-locked with the orbit (see Mamon 1987; 
Allen k Richstone 1988). 
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FIG. 3. Tidal stripping radius of clumps moving in different 

phase-locked circular orbits of radius Ri, for different mass fractions 
fji = m s /M c \. The initial clump size is 0.3 h' 1 Mpc. 

Figure 3 shows the tidal radius for such an orbit com- 
puted from King's (1962) 'instantaneous' criterion. Tidal 
stripping is negligible for fi > 0.05, except for very bound 
clumps (Ri < 0.6/i _1 Mpc). Note that during violent re- 
laxation, tides arising from transient density distributions 
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probably lead to gradients of the gravitational potential 
similar to those found in the central parts of relaxed clus- 
ters. Hence, clumps that do not survive to impulsive 
tides from the potential of a relaxed cluster would have 
hardly survived to the previous process of cluster merg- 
ing. Although the more massive clumps are more robust 
to tides, they are also more prone to dynamical friction, 
which causes their orbits to decay towards the cluster core 
where they disappear as substructures. 
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FIG . 4. Survival time (either to settle within R < R c or before 

being disrupted by tides) of detached cores vs. initial distance for 
radial (filled triangles) and circular orbits (filled squares), and for 
different mass fractions /i = m s /M c \. The dashed-dotted curves and 
horizontal dashed-dotted line denote the clump circular (upper) and 
radial (lower) unperturbed semi-orbital times and cluster crossing 
time (defined in § 3), respectively. The open triangle and square are 
the times spent within the cluster (R < Ri), before tidal disruption, 
by accreted groups, respectively on radial and non-radial orbits (see 
§3). 

Figure 4 shows the survival time (defined in § 2) for 
clumps of different mass. We adopt V% moving inwards and 
the same velocity dispersion for the cluster galaxies as in § 
3, yielding R l = 1.4 /i" 1 Mpc and M cl = 5.5 xlO 14 /^ 1 M Q 
for Slo = 1, but a larger value of <r g or Slo < 1 only cause 
slight changes in the results. When the clump survives 
the tides, (jjl > 0.07) the orbital decay time is shorter than 
an orbital period, and the clump falls directly into the 



cluster center. Nevertheless, clumps with 0.05 < [i < 0.2 
and Ri > 0.6 Mpc have a survival time longer than 
one cluster crossing time when their initial orbits are not 
too elongated (S > 0.3). Furthermore, only clumps with 
r s < 0.25 Mpc survive the cluster tides, which agrees 
with the upper limit to the sizes of detached cores derived 
from observation (SGS; GSS). 

5. DISCUSSION 

Small-scale substructure could be made of groups or 
small clusters that have been accreted since less than one 
cluster crossing time, and are about to be tidally dis- 
rupted. In this high rate of accretion would be 
required to reproduce the observed high fraction (roughly 
50%, § 1) of clusters with small-scale substructure. 

On the other hand, small-scale substructure may be 
caused by detached cores that have survived the violent 
relaxation after a merger of two similar sized clusters. Ac- 
cording to the survival times computed in § 4, the merger 
event could have taken place during the last 3 or 4 clus- 
ter crossing times for the longest-lived clumps, which are 
obtained by a fine-tuned balance between tides and dy- 
namical friction for /j, ~ 0.1-0.15 and initial circular or- 
bits. The probability that a detached core with /j, = 0.1 
has tgurv > ^cr is roughly 40%, if we assume that detached 
cores after violent relaxation are distributed with the same 
density profile as galaxies in the cluster (SGS), and with 
random initial orbit elongations. Therefore, the observed 
high fraction of detection of small-scale substructure may 
be reproduced, provided a high rate of cluster mergers. 

Nonetheless, the true origin of the observed small-scale 
substructure depends not only on the dynamical survival 
times, but also on the relative importance of the rates of 
accretion of small clusters versus the merging with clus- 
ters of similar mass, both of which depend on the power 
spectrum of the primordial density fluctuations and on 
Slo- In a forthcoming paper, we will compare the sur- 
vival timescales derived here with the observed frequency 
of small-scale substructure in apparently relaxed clusters 
to place important cosmological constraints. 

We thank the referee, Eliot Malumuth, for useful com- 
ments. This work was supported in part by CICYT 
projects PB 90-0448 and PB 92-0794. 
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